The two dimensional antiferromagnetic Heisenberg model in the presence of an 

external field 

M.S. Yang and K.H. Mutter 
Physics Department, University of Wuppertal, D-42097 Wuppertal, Germany 

(February 6, 2008) 



0\ 



m 
(N 

> 

0^ 
O 
O 



I 

o 
o 



We present numerical results on the zero temperature mag- 
netization curve and the static structure factors of the two 
dimensional antiferromagnetic Heisenberg model in the pres- 
ence of an external field. The impact of frustration is also 
studied. 

PACS number: 75.10 -b, 75.10.Jm, 75.40.Gb 



I. INTRODUCTION 

The discovery of high-Tc superconductorsEI has re- 
newed the interest in the two dimensional antiferromag- 
netic spin-i Heisenberg Hamiltonian : 



H= Six)Siy)+aJ2S{x)S{y). 



(1.1) 



<x,y> 



The ground state of (1.1) with nearest neighbors cou- 
pUngs <x,y> is considered to yield a good description 
of the antiferromagnetic properties in the undoped cop- 
per oxide planes. Doping destroys the autiferromagnetic 
order and this effect has been studiedo'cl intensively by 
adding to the nearest neighbor Hamiltonian a second 
term with next to nearest neighbor couplings [x, y\ in 
diagonal directions of strength a=J2/Ji- The behavior 
of the static structure factors Si{p = (pi,p2),a, ^) = 
A^exp(ipx)<Si{Q)Si{x)> - which are just the Fourier 
transform of the spin-spin correlators < Si{Q)Si{x) >, 
i=l,2,3 - changes with a : 

1. ) For < a < 0.3 one finds a strong divergence 

Si{p= {pi,p2),a,N oo) = ANm+\a) (1.2) 

at the momentum p = (tt, tt), indicating antiferromag- 
netic order with a nan vanishing staggered magnetization 
m+{a = 0) -0.308.0 

2. ) For a > 0.7 such a divergence is found a.t p — (0,7r) 
and p = (tt, 0), which is a signature for coUincar antifer- 
romagnetic order. 

It is unclear, what happens in between 0.3 < a < 0.7. 
One can imagine, that the singularity is less divergent 
and starts to move along some trajectory from p = (tt, tt) 
to p = (tTjO) and p = (0,7r). The reconstruction of this 
trajectory would demand a very precise determination of 



the static structure factor Si{p, a, N) with a high resolu- 
tion in the momentum plane. 

Moving singularities in static structure factors were ob- 
served and analyzed for the one dimensional antiferro- 
magnetifi-Heisenberg model in the presence of an external 
field B^B The relevant ground states, which enter in the 
computation of the static structure factors have total spin 
S3=S=M{B)N, where M^M{B) is given by the mag- 
netization curve. The following phenomena have been 
found in the one dimensional antiferromagnetic Heisen- 
berg model at fixed magnetization M : 

1) the longitudinal structure factor is almost constant 
for P3(M) = 7r(l - 2M) < p < tt : 

53(p,M) ^2(l-2Af) (1.3) 

and develops a cusp type singularity! ai p = p3{M). The 
critical exponent|473(Af ) at this singularity changes with 
the external fieldEl - i.e. with M. 

2) The transverse structure factor is almost constant for 
0<P <Pi{M) = 2'kM : 



Si{p,M) ~ 2M 



(1.4) 



and develops a "break" type singularity at the momen- 
tum pi{M). Moreover, it diverges at p = vr with a field- 
dependent critical exponent ?7i(M).EI 

The singularities in the static structure factors are 
accompanied with zero frequency excitations ("soft- 
modes") in the dynamical structure factors.&Ll The ef- 
fect of frustration on the position p3{AI), pi{M) and the 
type rj^LAI), rix{M) of the singularities has been studied 
as welhla Frustration does not change the position but 
drastically changes the type of the singularities. 

It is the purpose of this paper to investigate the singular- 
ities in the static structure factors of the two dimensional 
spin-i Heisenberg model in the presence of an external 
field. 

Our analysis is based on a numerical computation of the 
ground states at fixed magnetization M=S/N on the fol- 
lowing square lattices : 



7V = 4x4, iV = 6x6 
with periodic boundary conditions and 



(1.5) 



iV = 3x3+l, 7V = 5x5+l, 7V = 7x7+l (1.6) 
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with helical boundary conditionsB In the latter case the 
Hamiltonian can be considered to be one dimensional 



N 



H = J2 S{x){S{x + 1) + S{x + kj) 

x=l 

+ a I ^ S{x){S{x + k-l) + S{x + k + 1)) ) (1.7) 



with four types of couplings. Square lattices are realized 
for : 



iV = fc2 + 1, A: = 3,5,7 

= fc2 - 1, = 3,5,7 
as can be seen for fc=5, N—2Q in Fig. 1. 



(1.8) 



(1.9) 



The helical boundary conditions yield a quantization of 
the two dimensional momenta p — (pi,P2) : 
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(1.10) 



which differs from the usual quantization due to periodic 
boundary conditions. 

The outline of the paper is as follows. In section 2 we 
present the magnetization curves M(a, B) at zero tem- 
perature, as they depend on the external field B and 
the frustration parameter a. The field-dependence of the 
static structure factors at momentum p = (tt, tt) is dis- 
cussed in section 3. The momentum dependence of the 
static structure factors turns out to be smooth for the un- 
frustrated model, as is demonstrated in section 4. There- 
fore, we did not find any signature for field-dependent 
soft-modes in this case (a = 0.0). The situation changes 
if we frustrate the system sufficiently {a = 0.5). There we 
observe cusp-type singularities in the field-dependence of 
the static structure factors at certain momenta (section 
5). 



II. THE MAGNETIZATION CURVE 



Let us start with the ground state energies per site 
e{a, M=j^^ N) at fixed magnetization M=j^ and frus- 
tration a. The numerical results for the systems (1.5) 
and (1.6) scale with M except for the vicinity of M=0 
and M=^. Near saturation M ^ i, the ground state 
energy per site 



1 , a + 
e(a,M=-,7V) = — 
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and its first derivative 

6(a,M=i-^,7V) 



e(a,M=i,Ar) 



M - 



(2.1) 



-4, (2.2) 



with a < ^ are known for all system sizes N. 
Finite size effects appear in the difference ratios 



(M - i)-i(e(a, Af, N) - e(a, M = ^, iV)) 

^D{a,Z) 

= 4 - \D,{a)Z - \D2{a)Z^ - \D^{a)Z^ (2.3) 



They can be absorbed in the definition of an " optimized" 
scaling variable : 



Z =(l- 2AI) - 



c{M) 
N 



with 



c(M).2(iZ^)^ 



(2.4) 



(2.5) 



The choice for c(M) guarantees, that the expansion (2.3) 
satisfies the "initial" condition (2.2) for M=i — i.e. 
Z—0 independent of the exponent k. For k ~ 5, we 
achieve the best scaling of the numerical data as can be 
seen from Fig. 2. A fit of the numerical data to a Taylor 
expansion in Z yields the coefficients Dn{a), n = 1,2,3 
listed in table 1. We have repeated the fit with one ad- 
ditional term in the Taylor expansion (2.3). Such a fit 
yields only slight modifications of the first two coefficient. 

It is remarkable to note, that the first coefficient Di{a) 
decreases with a and vanishes at some value a=a* where 
0.49 < a* < 0.5. This has important consequence for the 
magnetization curve M=M(B, a) which is obtained from 
e(a, M, N = oo) by differentiation : 



de 
dM 



B{M). 



(2.6) 



For < a < a*, where -Di(a) > 0, we derive from (2.3) 
and (2.6) a linear behavior near saturation M — > i 



_ 1 A-B 
" 2 " iMa) 



(2.7) 



For a — a* - where Di{a* 
square root behavior : 



- however, we find a 



M 



1 



B 



2D2{a*) 



(2., 



A similar phenomenon has been observedo in the one di- 
mensional antiferromagnetic Heisenberg model with frus- 
tration a. There it turns out, that the magnetization 
curve has a square root behavior of the type (2.8) in the 
spin fiuid phase < a < ao with ~ 0.241 and a quar- 
tic root behavior in the vicinity of the transition point a, 
namely at a=0.25. 
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Let us go back to the discussion of the magnetization 
curves in the two dimensional model. The coefficients 
Dn{a), n = 1, 2, 3 in the polynomial (2.3) were fitted by 
the data points in the regime Z < ^ (i < A/ < i). 
The extrapolation of (2.3) to larger Z-values (smaller M 
values) is close to and away from the numerical data 
for a=0.0 and a=0.5, respectively. In other words : 
D{a — 0.0, Z) seems to be a smooth function in Z in 
contrast to D{a = 0.5, Z) which changes its slope sub- 
stantially somewhere between Z=0.3 and Z=QA. In or- 
der to determine the exact position Z, where the slope 
may be discontinuous, we would need numerical results 
on larger systems, which allow for a better resolution of 
the Z dependence of (2.3). A discontinuity in the slope 
of D{a = 0.5, Z) at Z = Zo (i.e. M = Mo) means, that 
the magnetization curve (2.6) has a plateau at height 
M = Mo. 

The existence of a plateau in the magnetizajtipn curve for 
a > 0.5 has been predicted by S. GluzmanllJ, who stud- 
ied the equivalent Bose gas problem and found a gap. 
In Figs. 3a, b we present the magnetization curves for 
a = 0.0 and a = 0.5 as they follow from an analy- 
sis of finite|-S|ystem results with the method of Bonner 
and Fisher O Finite size effects appear to be small for 
a = 0.0, as can be seen from a comparison of the nu- 
merical data for the small system sizes {N — 24, 26) with 
those for larger systems {N — 36,48,50) shown in the 
insets of Fig. 3a. 

For a = 0.5, however, we observe a sensitivity of the 
numerical data to the size of the system and to the 
boundary condition. The results for N = 24, 26 ob- 
tained with helical boundary conditions have a shoulder 
around A/ = i, which might indicate the emergence of 
a " plateau" . Unfortunately, we do not reach this shoul- 
der with our largest systems N = 48,50 shown in the 
upper insets of Fig. 3b. The lower inset of Fig. 3b 
contains the results for = 36 with periodic boundary 
conditions. Here, we observe a shoulder around A/=0.42 
which seems to be absent in the upper insets for the sys- 
tems with iV=48, 50. The variation of the magnetization 
curve with a on the 4_x 4 system has been computed by 
Lozovik and Notych.E3 These authors claim to see indi- 
cations for "plateaus" at various M-values for a ~ 0.538. 

III. THE FIELD-DEPENDENCE OF THE STATIC 
STRUCTURE FACTORS AT P=(7r,7r) 

In the one dimensional spin-i antiferromagnetic 
Heisenberg model, the field-dependence of the longitu- 
dinal and transverse structure factors at p=7r looks very 
differentcl : Finite size effects are small in the longitudinal 
case, the data points nicely scale with M (for M > 0) 
and develop a logarithmic singularity for M — > 0. Fi- 
nite size effects are very strong in the transverse case, 
the data points do not scale at all with M and indicate 
that the singularity persists at p—ir if the external field is 



switched on. The same behavior can be seen in the two 
dimensional case. In Fig. 4 we present the longitudinal 
structure factors S3{p — (tt, tt), Af, a. A) with a = 0.0 
and a = 0.5. All the data points follow one unique curve. 
Finite size effects turn out to be very small for larger val- 
ues of M. For decreasing values of M however, we find 
increasing finite-size effects. This is a signature for the 
emergence of a singularity in the limit Af ^ 0. Note also, 
that the longitudinal structure factor decreases with the 
frustration parameter a. 

The transverse structure factors Si{p — {tt,7t), AI,a, N) 
are shown in Fig. 5a and Fig. 5b for a^O.O and q;=0.5, 
respectively. The strong finite-size dependence for a=0.0 
indicates, that a weak magnetic field strengthens the sin- 
gularity at p={'k,tt) and thereby the antiferromagnetic 
order in the transverse structure factor. Note also the 
drastic change when we pass from the unfrustrated case 
(Fig. 5a) to the frustrated case (Fig. 5b). In the lat- 
ter case the finite-size dependence is less pronounced. 
Moreover, there appears a dip at Af=i in the data for 
A^ = 24,26 and at M=0A2 for A^ = 36 (inset Fig. 5b). 
In both cases, this is just the position, where we found a 
shoulder in the corresponding magnetization curves. 
The dip at Af=| has a good chance to survive in the 
thermodynamical limit; it is at least consistent with the 
numerical results on A^ = 24,26,36,48,50. The dip at 
M =0.42 in the data points seems to be a mere finite-size 
and/or boundary value effect. 

IV. THE MOMENTUM DEPENDENCE OF THE 
STATIC STRUCTURE FACTOR IN THE 
UNFRUSTRATED MODEL (q = 0) 

It was pointed out in the introduction that the static 
structure factors of the one dimensional spin-^ antiferro- 
magnetic Heisenberg model develop singularities at field- 
dependent momenta {pi{M)=2TrM, p3{M)=n{l - 2M)). 
We therefore address here the question, whether such sin- 
gularities can be found in the momentum dependence of 
the static structure factors for the two dimensional model 
as well. It will turn out, that the appearance of these sin- 
gularities crucially depends on the frustration parameter 
a. In this section we first treat the unfrustrated case 
a=0.0. 

The magnetization M=j^ is realized on three systems 
(1.5) (1.6) with N=24, 36, 48. 

The numerical values for the longitudinal structure fac- 
tor S'3(p= (pi,p2),M=^,a = 0, A), Ar=24, 36, 48 in the 
first Brillouin zone can be read off Fig. 6. These values 
appear to be constant along lines of constant ; 

a; = cospi -I- cosp2 (4-1) 

represented by the dashed curves in Fig. 6. The behav- 
ior can be seen directly along the line pi +P2=T^, x — 0. 
Therefore, we expect the longitudinal structure factor to 
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depend on the variable x only. In Fig. 7 we have plotted 
the longitudinal structure factor at M=j {N = 24,36), 
M=| {N = 24,48) and M=^ [N = 24,36,48) versus 
X. We observe a constant behavior for —2 < x < x^iM), 
and afterwards a decrease for X3{M) < x < 2. 2:3 (M) 
increases with M. 

The corresponding x-dependence of the transverse struc- 
ture factor is shown in Fig. 8. For x —2 {p ^ 
we see here the emergence of a singularity in accord with 
the strong finite-size dependence of the transverse struc- 
ture factor Si{p — {n,TT),M,a — 0,N) — 2M discussed 
in section 3. Moreover, this quantity is almost zero for 
X > xi{AI), where Xi{M) decreases with M. At M—0 
the longitudinal and transverse structure factor coincide. 
Their dependence on the variable (4.1) is shown in Fig. 
9 for N = 24, 26, 36. 

In the one dimensional case the longitudinal and trans- 
verse structure factors turned out to be constant for 
p > p3(M)=7r(l - 2M) and p < pi(M)=7r2M, respec- 
tively. However, there one finds cusp- type singularities at 
the "soft mode" momenta p=p3(M) andp=pi(M). The 
smooth momentum dependence of the static structure 
factors for the unfrustrated two dimensional Heisenberg 
model seems to exclude the existence of field-dependent 
soft-modes. 



V. HUNTING FOR SOFT-MODES 

In one and two dimensions, the antiferromagnetic 
Heisenberg models with nearest neighbor couplings differ 
substantially under various aspects : 

1. Fluctuations are stronger in one than in two dimen- 
sions i.e. the antiferromagnetic order increases with the 
dimension. 

2. There are field-dependent soft-modes in one dimen- 
sion; they seem to be absent in two dimensions. 

3. Near saturation (B ^ 4, M ^ i) the magnetiza- 
tion curves show a linear behavior in two dimensions but 
a square root behavior in one dimension. As was pointed 
out in section 2, this behavior changes, if the frustration 
parameter exceeds a critical value a*. 

If these features are correlated, we might hope to find 
soft-modes in the two dimensional antiferromagnetic 
Heisenberg model as well, provided we weaken the an- 
tiferromagnetic order by frustrating the system. In order 
to test this hypothesis, we need a reliable criteria for the 
existence or nonexistence of soft-modes. 

In the one dimensional easel, we looked for singulari- 
ties (breaks, cusps etc) in the momentum distribution 



of the structure factors at fixed magnetization M. As 
is demonstrated in Figs. 10a, b for the one dimensional 
case, pronounced structures are produced as well by the 
soft-modes in the M-dependence. The position of these 
singularities Mi{p) — ^ and M^lp) = define the 
soft-mode trajectories, which travel from p ^ (ferro- 
magnetic order) to p = tt (antiferromagnetic order) and 
vice versa. Keeping fixed the momentum, we have to 
meet the soft-mode at an appropriate value of the mag- 
netization. In this respect it is much easier to find the 
soft-modes in the one dimensional case. 
In case of the frustrated two dimensional Heisenberg 
model, the points p = (7r,7r) and p = (tt, 0), p = (0,7r) 
play a special role. Singularities in the structure factors 
at these points indicate antiferromagnetic and coUinear 
antiferromagnetic order, respectively. Therefore, it is 
rather plausible to assume, that the soft-mode trajec- 
tories connect these points. 

On the system with A^ = 5x5 — 1 = 24 sites, the 
momentum closest to p = (tt, tt) and p — (0, tt) are 
p — 7r(|,|) and p — 7r(i,|), respectively. Here, we 
have the best chance to observe soft-mode effects in the 
Af-distributions of the static structure factors. The sen- 
sitivity of the Af-distributions to a change of the frus- 
tration parameter can be seen in Figs. lla,b and 12. In 
the unfrustrated case a = 0.0 (open symbols), the M- 
distributions are smooth. Switching on the frustration 
parameter a (solid symbols), we observe the emergence 
of apeak at Af = i in Ssip = 7r(|, |), Af, a) (Fig. 11a) 
and at Af = i in S3{p = 7r(|, |), Af , a) (Fig. lib), re- 
spectively. The Af -distribution of the transverse struc- 
ture factor at p=(|, |)7r is shown in Fig. 12. The data 
for a = 0.5 show a break at Af=0.42, which reminds 
us to the structure observed in the one dimensional case 
shown in Fig. 10b. 

VI. DISCUSSION AND CONCLUSION 

In this note, we have studied the zero temperature 
properties of the two dimensional spin-^ antiferromag- 
netic Heisenberg model in the presence of an external 
field. Our results for the unfrustrated and the frustrated 
model can be summarized as follows. 

1. The unfrustrated model : Owing to small finite-size 
effects and a weak dependence on the boundary condi- 
tions, the magnetization curve and the momentum and 
field-dependence of the static structure factors can be 
determined with fairly good accuracy from the rather 
small systems (1.5), (1.6). This statement holds away 
from the "critical values" p—{TT,Tr), M—Q. At the mo- 
mentum p=(7r,7r) the longitudinal structure factor is fi- 
nite for Af > 0, whereas the transverse structure factor 
is divergent if the system size N goes to infinite. At 
fixed magnetization Af, the momentum dependence of 
the static structure factor is smooth and well described 
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by one variable, namely (4.1). There is no indication for 
the existence of singularities, which might give a hint to 
field-dependent soft-modes - i.e. zero frequency excita- 
tions. 



M. Karbach and G. Miiller for a critical reading of the 
manuscript and many useful comments. 



2. Frustration leads to new phenomena in the magne- 
tization curve and the static structure factors. The point 
a=0.b appears to be of special interest. Here, the mag- 
netization curve develops a square root singularity (2.8) 
near saturation (M For smaller values of M - at 

M=Mq{N) we find a shoulder in the magnetization curve 
and a dip in the field-dependence of the transverse struc- 
ture factor at p=(7r,7r). However, the position Mo{N) 
- where these phenomena appear - varies with the sys- 
tem size and the boundary conditions. We do not know, 
whether these phenomena are finite-size effects or will 
survive in the thermodynamical limit. 

We looked for singularities in the field-dependence of 
the static structure factors at fixed momenta. In- 
deed maxima emerge for momenta close to p=(7r,7r) 
and p=(0, 7r),p=(7r, 0) if the system is frustrated suffi- 
ciently. They might originate from a soft-mode with 
field-dependent momentum p{M) connecting p=(7r,7r) 
and p=(0, tt), p=(7r,0). However, our systems (1.5) (1.6) 
are much too small, in order to determine the field- 
dependence of the trajectory p{M). 

We have restricted our considerations to a limited range 
< a < ^ in the frustration parameter a. Going be- 
yond we meet level crossings. E.g. the ground 
states in the sectors with 3^ = ^ — ^ (called " one magnon 
states") have different momenta for a < ^ and a > i, 
respectively. Level crossings are "felt" by the Lanczos 
algorithm, which meets increasing problems to find the 
"true" ground state among two competing states. The 
problem can be solved, if these two states differ in their 
quantum numbers. One of the two competing states is 
filtered out, if the starting vector has the appropriate 
quantum numbers. Therefore, an extension of our anal- 
ysis to the strong frustration regime a > ^ demands for 
a careful analysis of the ground state quantum numbers 
in the sectors with a given total spin. We expect, that 
level crossings for a > ^ will generate interesting phe- 
nomena. E.g. they might be responsible for the plat&au 
in the magnetization curve predicted by S. GluzmanJ13 
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FIG. 1. The square lattice 5 x 5+1=26 with helical bound- 
ary condition 

FIG. 2. Scaling behavior of the difference ratio (2.3) in the 
strong field limit for the unfrustrated (a=0.0 open symbols) 
and the frustrated case {a=0.5 solid symbols). The curves 
represent polynomial fits in the scaling variable (2.4) 

FIG. 3. The magnetization curves computed on finite sys- 
tems with A'^=24,26,36,48 and 50 with the method of Bonner 
and Fisher [10] a) the unfrustrated case a=0.0 b) the frus- 
trated case with a=0.5 



FIG. 4. The field-dependence of the longitudinal structure 

factor S-dip = (hjTt), M,a, N) for a=0.0 (open symbols) and 
a=0.5 (solid symbols). The curves are shown to guide the eye 



FIG. 5. The field-dependence of the transverse structure 
factor Si{p= {n,n),M,a,N) for a) a=0.0, b) a=0.5 



FIG. 6. The momentum dependence of the longitudinal 
structure factor S3{p = {pi,p2),M = j2,a = 0.0, N) at fixed 
magnetization 

FIG. 7. Scaling of the longitudinal structure factor 
53 (p= {pi,p2), M,a = 0.0, N) in the variable a;=cospi+cosp2 
at fixed M=i,| and M=-^ 

FIG. 8. Same as Fig. 7 for the transverse structure factor. 
The curves are shown to guide the eye 



FIG. 9. Same as Fig. 7 for M=0.0 

FIG. 10. Field dependence of the one dimensional structure 
factors at fixed momenta p=|7r, |7r, ^tt, ^tt a) the longitu- 
dinal case b) the transverse case 

FIG. 11. Field dependence of the two dimensional longi- 
tudinal structure factors at fixed momenta a) p=7r(|, |) b) 
p=7r(i, |) open symbols a=0.0 and solid symbols q;=0.5 

FIG. 12. Field dependence of the two dimensional trans- 
verse structure factor at fixed momentum p=7r(|, |) 



TABLE I. The coefficients £>„(«) of the expansion (2.3) 

a 0.3 0.4 0.45 0.49 0.495 0.5 

Dija) 1.653 1.167 0.802 0.478 0.063 -0.003 -0.074 

D-ija) 1.374 1.086 1.162 1.763 3.115 3.376 3.663 

Dala) 0.918 2.928 3.077 2.317 -0.012 -0.482 -1.003 
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